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1. Introduction 

VER J0521+211, spatially associated with the radio and X-ray source RGB J0521.8-1-2112, 
was discovered in 2009 as TeV emitter by VERITAS [1] and is now unambiguously classified 
as a BE Eac-type blazar [2]. As usually observed in blazars, VER J0521-I-211 is highly vari¬ 
able at all wavelengths. With an average integral flux above 200 GeV of about >0.1 Crab^, 
VER J0521-I-211 is among the brightest known TeV blazars. Eocated at a redshift of z = 0.108 [4], 
this flux corresponds to a gamma-ray luminosity of about 3 x lO^'^ergs^^ and is thus significantly 
brighter than the classical northern high-synchrotron-frequency-peaked BE Eacs (HBEs) such as 
Mrk 421, Mrk 501 and lES 1959-1-650. VER J0521-I-211 itself shows a synchrotron peak in the op¬ 
tical regime, suggesting a classification as intermediate-frequency-peaked BE Eac (IBE). However, 
during the TeV high state in November 2009 [5], VER J0521-I-211 showed more HBE-like prop¬ 
erties, especially given the observed spectral hardening in the X-ray regime with increased flux. 
The observed variability and its bright TeV flux motivated further multiwavelength observations 
of VER 1052 1-1-2 11 in order to extend the detailed time-resolved spectral modeling available for 
nearby HBEs [6, 7, 8] to a more luminous blazar which are reported here. 

2. VERITAS 

VERITAS is an array of four 12 m-diameter imaging atmospheric Cherenkov telescopes (lACTs), 
located in southern Arizona, USA. It is sensitive to gamma-ray energies between 85 GeV and 
>30 TeV and is capable to detect a source with 0.01 Crab flux in under ~ 25 hours of observa¬ 
tions. VERITAS observations of VER 10521-1-211 reported here were taken in two seasons: the 
2012 season comprising data taken between 2012 Nov 11 and 2013 Eeb 14; and the 2013 season 
(2013 Oct 15 to 2014 Eeb 01). All data were taken in wobble mode for which the source is ob¬ 
served with an offset of 0.5° with respect to the center of the instrument’s field of view to allow for 
simultaneous background measurements [9]. After quality selection, the dataset comprised 5.0h 
live time in the 2012 season and 16.5 h in the 2013 season. VERITAS data analysis followed the 
procedure outlined in [10]. Cherenkov light from air showers initiated by gamma rays and cosmic 
rays triggers the readout of PMT signals, which are then calibrated and used to reconstruct an im¬ 
age of the shower in the focal plane. Individual telescope images are parameterized by moment 
analysis, and geometrical reconstruction is used to estimate the arrival direction, 0, of the primary 
particle with respect to the source location. The signal was extracted using a reflected region back¬ 
ground model with an ON region of 0.17° radius centered on the position of VER 10521 -i-211 (R. A. 

= 05’’21™45®, Deck = -|-21°12'51."4, 12000). The strength of the gamma-ray signal is evaluated 
using Eq. 17 in [11] and was found to be 19.0 a in 2012 and 60.3 a in 2013. The average flux levels 
in the two seasons are 7Ao^2TeV ~ 2.4 ±0.2 andF^Q^ 2 TeV ~ 4.2 ±0.1 in units of 10^^' cm^^s^^ Eor 
comparison, the time averaged flux in [2] was 1.9 ±0.1 in the same units (as shown in Eigure 1). 

3. MAGIC 

MAGIC is a stereoscopic system of two 17 m-diameter lACTs, located at Ea Palma, Canary 

*The gamma-ray flux of the Crab Nebula used here is F{E > 200GeV) = 2.1 x 10^*°cm^^s^* [3]. 
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Figure 1: Evolution of the TeV flux. X-ray flux, and X-ray hardness ratio as a function of time during the 
2012 {left) and 2013 seasons {right). It should be noted that the MAGIC and VERITAS observations in 2013 
were not taken simultanously. 

Islands. It is sensitive to gamma-ray energies above ~ 50 GeV (in normal trigger mode) and under¬ 
went a series of upgrades during summer 2011 and 2012 [12]. In stereoscopic mode, the integral 
sensitivity for sources with Crab Nebula-like spectra above 220 GeV is (0.66 ±0.03)% of Crab 
Nebula flux in 50 h of observations (more details about the MAGIC performance are given in [13]). 
MAGIC observations of VER J0521-I-211 were taken between 2013 Oct 15 and Dec 2 and were 
triggered by the high state in optical and high-energy gamma rays reported by KVA (E. Eindfors, 
private communication) and Fenni-LKT [14], respectively. Observations covering a zenith angle 
range from 7° to 26° were performed in wobble mode with an offset of 0.4°. After quality selec¬ 
tion and dead time correction, the total effective observation time was ~ 4.5 hours. The data were 
analyzed using the MAGIC analysis and reconstruction software package [15], which has been 
adapted to stereoscopic observations [16]. The signal was extracted from an ON region of 0.11° 
radius centered on the position of VER J0521-I-211. The analysis of the overall data set yielded a 
strong signal of 30.5 a significance [11], while the significances of the signal from the individual 
nights are 18.6 a, 9.7 a, 21.9 a and 6.2 a for Oct 15, 16, Nov 29 and Dec 2, respectively. Given 
the significance of the signal, night-wise spectra could be obtained above the energy threshold of 
65 GeV for the observations in October and November and are shown in Eigure 2. 

4. X-ray 

The Swift X-ray telescope (XRT [17]) monitored VER J0521-I-211 in the X-ray band (0.3- 
10 keV) for most of the nights of the campaign. Data were processed using the most recent versions 
of the standard Swift tools: Swift Software version 3.9, ETOOES version 6.12 and XSPEC version 
12.7.1. Eight curves are generated using xrtgrblc version 1.6. More details on the monitoring 
program and the analysis procedure can be found in [18]. 

5. Flux and spectral variability 

The TeV fluxes measured in 2012 and 2013 are higher than those reported during the discovery 
of the source [2], with the average flux in 2013 being a factor of 2 brighter than that of 2009. The 
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Figure 2: Spectral energy distribution of VER J0521+211 as obtained by MAGIC during three single nights 
in 2013. 
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Figure 3: Comparison of the spectral energy distributions measured with VERITAS at different flux states. 
Spectra have been obtained by stacking observations with nightly fluxes F>o,2TeV <3x10^** cm^^ s^ * (low 
flux), F>o,2TeV > 4 X 10^" cm^^ s^* (high flux), and a medium flux bin defined between these two values. 
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nightly flux (Figure 1) shows significant variability down to timescales of 1 day. 

To test for spectral variability, we have produced gamma-ray spectra for individual observing 
nights (Figure 2) and by stacking observations at different flux levels (Figure 3). At low and 
medium flux levels we find gamma-ray spectra largely compatible with a power-law with index 
F = 3.44±0.20 that was reported in previous observations [2]. Co-adding the spectra from the three 
nights with higher fluxes we measure F = 3.10 ± 0.07, constituting evidence for a mild hardening 
of the gamma-ray spectrum when the emission is brighter. In addition, the lower threshold of the 
MAGIC observations (Figure 2) helps to reveal some evidence of spectral curvature, suggesting a 
peak of the inverse-Compton emission component at energies of ~ 100 GeV. 

Similar levels of 1-day-scale flux variability can be seen in the X-ray band (Figure 1). The 
evolution of the hardness ratio in the X-ray band (Figure 4) suggests only small changes of the 
X-ray spectral shape at different flux levels. Data from the nights with simultaneous observations 
by VERITAS and Swift-XRI (Figure 5) show a clear correlation between the measured fluxes in 
both bands. 



Figure 4: X-ray hardness ratio as a function of X-ray flux measured by Swift-XKT. The hardness ratio is 
defined as the rate in the 2-10 keV band divided by the rate in the 0.3-2 keV band. 


6. Discussion 

The data presented in these proceedings constitute the first extensive multi-wavelength cam¬ 
paign on VER J0521-I-211. Significant changes in the X-ray and TeV fluxes from the source were 
seen in timescales of ~ 1 day, with an amplitude of approximately a factor of 5 in both bands. 
Spectral variability in the TeV band is first reported here, displaying a mild hardening of the spec¬ 
trum at higher flux levels. The spectral variability is however much less significant than typically 
found in HBLs. 

The observed correlation between the fluxes in the X-ray and TeV bands suggests a connection 
between the synchrotron and the high-energy spectral components of the emission, as expected 
if the high-energy emission is dominated by synchrotron self-Compton. However, a contribution 
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Figure 5: X-ray flux as a function of the TeV flux for nights with simultaneous observations in both bands. 


from inverse-Compton emission from electrons and positrons interacting with other radiation fields 
external to the jet cannot be excluded. 

Compared to similar studies conducted on TeV HBLs (Mrk 421, Mrk 501, lES 1959-1-650), the 
level of spectral variability in the TeV regime is small given the high-significance of fhe defecfions 
and fhe TeV flux changing by a factor of ~ 5. This small level of specfral variabilify and a softer 
specfral index seem characferisfic of IBLs, where elecfrons and posifrons wifh energies beyond fhe 
TeV scale will be less dominanf fhan in HBLs given fhe higher rate of inverse-Compfon cooling. 
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